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ABSTRACT 


The  Apollo  Lunar  Program  utilized  efficient,  i.e..  Earth-return,  transearth 
trajectories  which  employed  parking  orbits  in  order  to  minimize  energy  requirements. 
This  thesis  concentrates  on  a  different  type  of  transearth  trajectory.  These  are  direct- 
ascent,  hyperbolic  trajectories  which  omit  the  parking  orbits  in  order  to  achieve  short 
flight  times  to  and  from  a  future  lunar  base.  The  object  of  the  thesis  is  the  development 
of  a  three-dimensional  transearth  trajectory  model  and  associated  computer  program  for 
exploring  trade-offs  between  fiight-time  and  energy,  given  various  mission  constraints. 
The  program  also  targets  the  Moon  with  a  hyperbolic  trajectory,  which  can  with  a  time- 
reversed  trajectory;  be  used  for  targeting  Earth  impact  points.  The  first-order  model  is 
based  on  an  Earth-centered  conic  and  a  massless  spherical  Moon,  using  MathCAD 
version  3.0.  This  model  is  intended  as  the  basis  for  future  patched-conic  formulation  for 
the  design  of  fast  Earth-return  trajectories.  Applications  include  placing  nuclear- 
deterrent  arsenals  on  the  Moon,  various  space  support  related  activities  and  finally 
protection  against  Earth-threatening  asteroids  and  comets  using  lunar  bases. 


iii 


TABLE  OF  CONTENTS 


I.  INTRODUCTION  . 1 

A.  BACKGROUND  .  1 

B.  OBJECTIVES  AND  LIMITATIONS .  2 

II.  ANALYSIS  APPROACH .  3 

A.  TRANSLUNAR  TRAJECTORY  DESIGN .  4 

1.  Translunar  Trajectory  Types  .  5 

2.  Non-Inertial  Coordinate  Frames  .  5 

B.  SOLUTION  FOR  POINT-MOON  INTERCEPT .  7 

1.  Three  Dimensional  Point-Moon  Model  .  8 

2.  Three  Dimensional  Considerations .  8 

C.  SOLUTION  FOR  MASSLESS  SPHERICAL  MOON  INTERCEPT  .  .  12 

1.  Iterative  Technique  for  Lunar  Surface  Impact  .  12 

2.  Selenographic  Coordinates  .  14 

a.  Euler  angles .  14 

b.  Lunar  Impact .  14 

D.  PATCHED  CONIC  APPROXIMATION .  17 

III.  APPLICATIONS .  19 


iv 


A.  STRATEGIC  PEACE  INITIATIVE  CONCEPT  .  19 

1.  Background  Information .  19 

2.  SPI  Architecture . 21 

B.  SPI  TECHNICAL  FEASIBILITY  . 22 

C  SPACE  SUPPORT  ACTIVITIES  . 24 

D.  SUMMARY . 24 

IV.  CONCLUSIONS  AND  RECOMMENDATIONS . 26 

A.  FINDINGS . 26 

B.  RECOMMENDATIONS . 26 

APPENDIX  A.  PATCHED-CONIC  BASED  ON  INERTIAL  FRAMES  . 27 

APPENDIX  B.  LUNAR  TRAJECTORY  PROGRAM  . 28 

LIST  OF  REFERENCES . 55 

INITIAL  DISTRIBUTION  LIST . 56 


V 


TABLE  OF  SYMBOLS  AND  ABBREVIATIONS 


This  listing  of  symbols  and  abbreviations  is  to  aid  the  reader  in  clarifying  the 
multitude  of  symbols  and  abbreviations  used  in  this  thesis.  They  are  in  alphabetical  order. 

“e  orb  =  input  value  for  posigrade  or  retrograde  orbit 

ag  :  spacecraft  right  ascension  at  Earth  launch 

Om  :  spacecraft  right  ascension  at  Moon  intercept 

ckm,  :  represents  cosine  of 

ofMj  :  represents  sine  of 

<*EM  ■  represents  the  difference  between  and 

«EMi  •  represents  the  cosine  of 

«EM2  •  represents  the  sine  of 

Oeo  :  represents  the  difference  between  the  translunar  trajectory  longitude-of- 
ascending-node  and  spacecraft  right  ascension  at  Earth  launch 
«Ein  :  represents  the  cosine  of  agn 
“En2  •  represents  the  sine  of 

a,  :  represents  an  intermediate  value  to  calculate  the  spacecraft  right 

ascension  on  its  trajectory 

Oq  :  represents  spacecraft  right  ascension  on  its  trajectory 

a  ;  represents  the  selenographic  longitude  for  lunar  surface  impact 

/Sg  :  represents  flight  path  azimuth  at  launch 
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Pei 

:  represents 

Pe2 

:  represents 

Pm 

:  represents 

Pmi 

:  represents 

Pm2 

:  represents 

Po 

:  represents 

Px 

:  represents 

P2 

:  represents 

:  represents 

:  represents 

Ce2 

:  represents 

:  represents 

Co 

:  represents 

Ye 

:  represents 

Yei 

:  represents 

Ye2 

:  represents 

Ym 

:  represents 

Ymi 

:  represents 

YM2 

:  represents 

Yo 

:  represents 

Yi 

:  represents 

Y2 

:  represents 

D 

:  represents 

cosine  of 
sine  of/S,. 

flight  path  azimuth  at  moon  intercept 
cosine  of 
sine  of 

flight  path  azimuth  of  the  spacecraft  on  its  trajectory 
cosine  of  )3o 
sine  of  /S,, 

the  sweep  angle  at  launch 
the  cosine  of 
the  sine  of 

the  sweep  angle  at  moon  intercept 

the  sweep  angle  of  the  spacecraft  on  its  trajectoiy 

flight  path  angle  at  launch 

the  cosine  of  ye 

the  sine  of  Ye 

flight  path  angle  at  moon  intercept 
the  cosine  of  Ym 
the  sine  of  Ym 

the  flight  path  angle  of  the  spacecraft  on  its  trajectory 

the  cosine  of  yo 

the  sine  of  Yo 

the  deflection  angle 
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5e 

Sm 

So 

5 

Aa 

A«m 

AaM2 

Atto 

AO) 

Aa2 

^pO 

Ah, 


eccen 


eccen„ 


energy„i„ 


:  represents  the  directrix-to-focus  or  axial  offset  of  the  asymptotic  from 
the  focus  ;  it  is  a  variable  valid  for  eccentricities  greater  than  one 
:  represents  the  distance  between  the  spacecraft  and  the  lunar  center  at 
time  to 

;  represents  the  spacecraft  declination  at  earth  launch  point 
;  represents  spacecraft  declination  at  moon  intercept 
:  represents  the  spacecraft  declination  on  its  trajectory 
:  represents  the  selenographic  latitude  for  surface  impact  on  the  moon 
;  see  Figure  8 

:  represents  intermediate  value  to  calculate  spacecraft  right  ascension  on 
its  trajectory 

:  represents  the  cosine  of  Aa^ 

:  represents  the  sine  of  Ao-m 

;  represents  an  intermediate  value  for  the  calculation  of  Oq 
:  represents  the  cosine  of  Aao 
:  represents  the  sine  of  Aao 

:  represents  cosine  of  ;  see  iterative  process  on  page  12 
;  represents  the  difference  between  point-moon  intercept  flight  time 
relative  to  perigee  passage  and  the  spacecraft  flight  time  on  its 
trajectory 

:  represents  eccentricity  as  an  input  value 
:  represents  the  minimum  eccentricity  input  value  allowed 
:  represents  the  minimum  specific  energy  required  for  launch 
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Eh,  :  represents  the  hyperbolic  eccentric  anomaly  for  earth  launch 

Eh2  :  represents  the  hyperbolic  eccentric  anomaly  for  point-moon  intercept 

Eh3  :  represents  the  hyperbolic  eccentric  anomaly  for  the  spacecraft  on  its 

trajector>' 

E,hi  :  represents  the  hyperbolic  sine  of  Eh, 

E.h2  ■  represents  the  hyperbolic  sine  of  Eh2 

E,h3  :  represents  the  hyperbolic  sine  of  Ehj 

fy,  :  represents  the  input  value  for  the  true  anomaly  of  the  spacecraft  at 

moon  intercept 

fg  :  represents  the  true  anomaly  of  the  spacecraft  at  earth  launch 

;  represents  the  difference  of  the  true  anomaly  of  the  spacecraft  at  moon 
intercept  and  at  earth  launch 

MAX  •  represents  a  constraint  value  on  fjv,  if  required  to  be  used 
fM.MAxi  :  represents  the  cosine  of  f^j^^ 
fM_MAX2  represents  the  sine  of  I'm  max 
fh,,  :  represents  the  cosine  of  f^ 

fM2  :  represents  the  sine  of  f^ 

fo  :  represents  the  true  anomaly  of  the  spacecraft  on  its  trajectory 

^MO  :  represents  '  he  cosine  of  f^ 

0  :  represents  the  first  Euler  angle 

he  :  represents  the  geocentric  altitude  of  launch  point 

i^,  :  represents  the  inclination  of  the  moon’s  orbit  plane 

i  :  represents  the  translunar  trajectory  inclination 
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i,  :  represents  the  cosine  of  i 

i,  :  represents  the  sine  of  i 

kg  :  represents  the  earth’s  gravitational  constant 

kv,  :  represents  the  moon’s  gravitational  constant 

rg  :  represents  the  earth’s  mean  radius 

r^  ;  represents  the  moon’s  mean  radius 

Re,  :  represents  the  spacecraft  radial  distance  at  earth  launch 

^Ml  :  represents  the  spacecraft  radial  distance  at  moon  intercept 

r  :  represents  the  radial  distance  ratio 

:  represents  apogee  distance  :  valid  only  for  eccentricities  less  than  1 
Rp  :  represents  the  perigee  distance 

Rq  :  represents  the  input  value  for  the  position  of  the  spacecraft  on  its 

trajectory 

Rm-  :  represents  the  radial  velocity  at  point-moon  intercept 
R,  :  represents  an  intermediate  value  to  calculate  fp 

Pem  •  represents  the  Earth-Moon  mean  distance 

sigOjEM  •  represents  an  intermediate  value  calculated  to  insure  correct  use  of  sign 
when  calculating  aEM 

SfEM  •  represents  the  input  value  of  + 1  or  -  1  to  insure  correct  sign  usage  in 
calculating  Oem 

Tp  :  represents  a  circular  orbit  period  for  radius  Rp 

tg  :  represents  Earth  launch  flight  time  relative  to  perigee  passage 

tg,  :  represents  an  intermediate  value  used  to  calculate  tg 
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:  represents  the  point-moon  intercept  flight  time  relative  to  perigee 
passage 

^M)  •  represents  an  intermediate  value  used  to  calculate  t^ 

•  represents  the  point-moon  translunar  trajectory  flight  time 

Imo  :  represents  an  intermediate  value  used  to  calculate  t,, 

to  :  represents  the  flight  time  of  the  spacecraft  on  its  trajectory 

Vp  ;  represents  the  perigee  velocity 

Vg  :  represents  the  spacecraft  velocity  at  launch 

:  represents  the  spacecraft  velocity  at  moon  intercept 
Vo  :  represents  the  spacecraft  velocity  on  its  trajectory 

•  represents  the  input  value  of  the  lunar  sweep  angle  at  moon  intercept 

<f>^^  :  represents  the  cosine  of 

<t>M2  •  represents  the  sine  of 

0^,0  •  represents  the  difference  between  the  lunar  sweep  angle  at  moon 

intercept  and  the  moon’s  mean  orbital  rotation  rate  multiplied  by 
4)  :  represents  the  first  Euler  angle 

Xm  :  represents  the  X-coordinate  for  the  rectangular  coordinates  of  the  moon 

at  point-moon  intercept 

Xm-  ■<  represents  the  X-component  for  calculating  radial  velocity  at  point-moon 
intercept 

Xq  :  represents  the  X-coordinate  for  the  rectangular  coordinates  of  the 

spacecraft  on  its  trajectory 

Xmo  :  represents  the  X-coordinate  for  the  center  of  the  moon 
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X  :  represents  the  difference  between  Xm  and  X^o 

^lun  •  represents  the  X-coordinate  for  the  rectangular  selenographic 
coordinates 

:  represents  the  Y-coordinate  for  the  rectangular  coordinates  of  the  moon 
at  point-moon  intercept 

Y^.  :  represents  the  Y-component  for  calculating  radial  velocity  at  point-moon 

intercept 

Yq  :  represents  the  Y-coordinate  for  the  rectangular  coordinates  of  the 

spacecraft  on  its  trajectory 

Ymo  •  represents  the  Y-coordinate  for  the  center  of  the  moon 

Y  :  represents  the  difference  between  Y^  and  Y^o 

Y,u„  ;  represents  the  Y-coordinate  for  the  rectangular  selenographic 
coordinates 

Tm  :  represents  the  third  Euler  angle 

Zm  :  represents  the  Z-coordinate  for  the  rectangular  coordinates  of  the  moon 
at  point-moon  intercept 

Zm  :  represents  the  Z-component  for  calculating  radial  velocity  at  point-moon 
intercept 

Zmo  represents  the  Z-coordinate  for  the  center  of  the  moon 

Z  :  represents  the  difference  between  Z„  and  Z^o 

Ziun  :  represents  the  Z-coordinate  for  the  rectangular  selenographic 
coordinates 

n  ;  represents  the  translunar  trajectory  longitude-of-the-ascending-node 
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fc>E  :  represents  the  Earth’s  angular  rotational  rate 

(Om  :  represents  the  Moon’s  mean  orbital  rotation  rate 

G)  :  translunar  trajectory  argument-of-perigee 
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I.  INTRODUCTION 


This  chapter  provides  a  general  description  of  this  thesis  and  is  divided  into  two 
sections.  The  first  section  presents  background  information  pertaining  to  translunar 
trajectory  design.  The  second  section  delineates  the  objectives  and  limitations  of  this 
project. 

A.  BACKGROUND 

While  the  impending  mission  of  establishing  a  permanent  lunar  base  stands  in  our 
future,  it  is  appropriate  to  study  this  type  of  mission  by  investigating  various  translunar 
trajectories.  The  trade-off  between  (light  time  and  energy  required  must  be  determined 
prior  to  considering  the  type  of  translunar  trajectory  to  design.  The  energy  required  for 
the  hyperbolic  trajectory  can  then  be  compared  to  that  of  the  Apollo  mission  to 
determine  the  increase  in  energy  and  thereby  cost  increases. 

The  Apollo  Lunar  Program  of  the  1960s  utilized  energy-efficient  translunar 
trajectories.  This  thesis  will  analyze  a  different  class  of  trajectories  named  fast  Earth- 
return  trajectories.  These  trajectories  utilize  a  direct-ascent  from  the  lunar  surface  and 
omit  the  lunar  parking  orbit  used  in  the  Apollo  missions.  These  direct-ascent  trajectories 
differ  considerably  from  translunar  trajectories  previously  considered  and  discussed  in  this 
thesis.  (Wadsworth,  1991) 
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B.  OBJECTIVES  AND  LIMITATIONS 

The  objective  of  this  thesis  is  to  develop  a  trajectory  model  and  associated 
computer  program  for  systems  studies  of  fast  Earth-return  trajectories  from  a  lunar  base. 
The  scope  of  the  project  involves  a  step  process  which  evolves  in  complexity,  as  the 
design  moves  from  a  massless  point-Moon,  to  a  massless  spherical  Moon  intercept 
solution.  Finally,  the  groundwork  is  laid  for  a  patched-conic  model  which  culminates  in 
a  solution  for  intercept  on  a  Moon  with  finite  lunar  mass.  A  novel  feature  is  the  use  of 
a  time-reversed  trajectory  for  targeting  Earth  impact  points.  In  other  words,  the  lunar 
impact  point  serves  as  the  lunar  launch  point  for  a  time-reversed  trajectory.  The 
complete  patched  conic  approximation,  as  discussed  in  Chapter  II,  is  not  developed  in  the 
program  due  to  scope  of  work  and  time  constraints. 
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II.  ANALYSIS  APPROACH 


This  chapter  sets  the  foundation  of  the  thesis  by  delineating  a  flow  chart  of  tasks 
that  are  required  to  be  accomplished  to  reach  the  objective  of  the  thesis.  A  flow  chart 
represented  in  Figure  1,  summarizes  the  work  to  be  accomplished.  The  last  block  of  the 
diagram  represents  further  research  work. 


Figure  l.  Translunar  Trajectory  Design  Flowchart. 

Precision  translunar/transearth  trajectory  design  requires  time-consuming,  step-by- 
step  numerical  solution  of  the  exact  differential  equations  of  motion,  including  lunar  and 
solar  ephemerides.  For  preliminary  design,  the  computation  time  can  be  reduced  a 


hundred-fold  or  more  by  employing  a  closed-form  "patched-conic"  approximation  to  the 
precision  trajectory.  Since  the  patched-conic  is  a  composite  of  i;he  two  conics,  it  can  be 
formulated  explicitly  in  terms  of  elementary  functions  oi  the  key  parameters.  This 
provides  better  insight  into  design  trade-offs  than  numerical  computation.  The  patched 
conic  approximates  the  gravitational  effects  of  the  Earth-Moon  system.  It  ignores  the 
small  perturbation  due  to  the  Sun. 

Two  approaches  were  considered  for  developing  the  patched-conic  model.  Both 
approaches  employ  Earth-centered  and  Moon-centered  coordinate  frames  with  respect 
to  which  the  conic  trajectories  are  formulated.  Although  the  conic  trajectories  are  planar 
with  respect  to  the  Earth-centered  and  Moon-centered  Frames,  they  form  three- 
dimensional  curves  as  viewed  in  an  inertial  barycentric  frame.  In  the  first  approach,  the 
coordinate  frames  are  non-inertial,  being  attached  to  the  Earth  and  Moon  which  travel 
in  (idealized)  circular  orbits  (with  respect  to  the  Earth-Moon  inertial  frame).  The  second 
approach  employs  non-rotating,  inertial  Earth-centered  and  Moon-centered  coordinate 
frames.  These  frames  are  attached,  respectively,  to  a  fictitious  Earth  and  fictitious  Moon 
which  travel  at  constant  velocity  with  respect  to  the  barycentric  frame.  The  implications 
of  the  second  approach  are  briefly  stated  in  Appendix  A.  The  advantage  of  the  first 
approach,  pursued  in  this  thesis  research,  is  described  in  section  A.  (Wadsworth,  1991) 

A.  TRANSLUNAR  TRAJECTORY  DESIGN 

The  first  consideration  given  to  designing  a  translunar  trajectory  is  to  decide  on 
what  translunar  trajectory  types  will  be  used.  Since  this  thesis  concentrates  on  fast  Earth- 
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return  trajectories,  eccentricities  greater  than  one  will  be  utilized.  The  reason  for  using 
non-inertial  coordinate  frames  is  discussed  in  section  A.2, 


1.  Transhinar  Trajectory  Types 

Figure  2  illustrates  all  possible  trajectory  types  that  can  be  considered  for 
Earth-Moon  trajectory  design.  The  radius  of  the  Earth  (Rn),  radius  of  the  Moon  (Rm), 
spacecraft  true  anomaly  at  launch  (f^),  and  spacecraft  true  anomaly  at  moon  intercept 
(f^,)  are  projected  in  Figure  2. 

Since  only  trajectories  with  eccentricities  greater  than  1.0  are  being 
considered,  only  type  I  and  III  hyperbolic  trajectories  will  be  used  in  thesis  research  (see 
Table  I). 

2.  Non<lnertial  Coordinate  Frames 

The  Earth-centered  and  Moon-centered  coordinate  frames  are  non-inertial 
frames,  since  they  are  in  a  state  of  gravitational  free-fall  about  the  Earth-Moon 
barycenter  (center  of  gravity).  As  a  consequence,  the  composite  patched-conic  provides 
a  better  approximation  to  the  restricted,  three-body  spacecraft  motion  than  would  be 
suggested  by  simply  joining  two  conics. 

Within  the  Moon’s  region-of-influence,  the  terrestrial  perturbing  acceleration 
of  the  spacecraft  is  nearly  the  same  as  the  centripetal  acceleration  of  the  Moon-centered 
coordinate  frame.  In  this  frame,  these  accelerations  nearly  cancel;  this  can  be  deduced 
from  the  facts  that  the  distances  of  the  spacecraft  to  the  Earth  and  the  Moon  to  the 
barycenter  are  comparable  and  the  mass  of  the  Earth  is  comparable  to  the  equivalent 
mass  at  the  barycenter.  For  the  Earth-centered  phase,  there  is  a  similar  benefit  over 
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Figure  2.  Equatorial  plane  projection. 


TRANSLUNAR  TRAJECTORY  TYPES: 

TYPE  I  poaigTsde  :  arc  ab 
TYPE  I  retrograde  :  arc  dc 
TYPE  II  poaigrade  :  arc  abc  (  e  <  1  ) 

TYPE  II  retrograde  ;  arc  deb  (  e  <  1  ) 

TYPE  in  poaigrade  ;  arc  dab 
TYPE  in  retrograde  :  arc  adc 
TYPE  rV  poaigrade  :  arc  dabc  (  e  <  1  ) 

TYPE  rV  retograde  ;  arc  adeb  (  a  <  1  ) 


Table  I.  Trajectory  Types. 
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most  of  the  trajectory,  especially  for  parabolic  and  hyperbolic  trajectories.  This 
cancellation  does  not  occur  in  the  case  of  patched-conics  based  on  inertial  frames  as 
described  in  Appendix  A  (Wadsworth,  1965). 

The  Earth-centered  translunar  conic  trajectory  may  be  either  elliptic  or 
hyperbolic.  The  latter  corresponds  to  shorter  flight  times  and  is  the  only  case  that  needs 
consideration  for  fast  trajectories.  For  the  same  reason,  the  moon-centered  conic 
approach  trajectory  is  also  modeled  as  hyperbolic  in  this  study.  Both  the  Earth-centered 
and  Moon-centered  frames  are  non-rotating.  The  motions  of  the  Earth  and  Moon  about 
their  common  barycenter  are  approximated  by  coplanar  circular  orbits.  In  general,  the 
plane  of  these  circular  orbits  does  not  coincide  with  the  three  dimensional  trajectory  of 
the  spacecraft.  (Wadsworth,  1991) 

B.  SOLUTION  FOR  POINT-MOON  INTERCEPT 

The  least  complicated  form  of  a  first  order  solution  is  for  a  point-moon  intercept. 
Figure  3  illustrates  the  relevant  parameters  of  the  problem.  Appendix  B,  utilizing 
MathCAD  3.0,  numerically  solves  the  first  order  solution  of  point-Moon  and  massless 
spherical  Moon  intercept.  In  Figure  3  the  following  assumptions  are  made: 

•  The  Earth-centered  cartesian  coordinate  frame  is  non-rotating  (indicated  by 
X,Y,and  Z). 

•  The  Moon-centered  cartesian  coordinate  frame  is  non-rotating  (indicated  by  x,y,and 
z  with  the  x-axis  parallel  to  the  X-axis,  and  the  z-axis  tilted  at  angle  for  the  Z- 
axis.) 
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The  initial  point-Moon  intercept  solution  represents  a  simple  three-dimensional 
model;  three-dimensional  considerations  arc  incorporated  in  the  first-order  solution  solved 
in  Appendix  B.  These  considerations  will  be  discussed  in  section  B.2. 

1.  Three  Dimensional  Point-Moon  Model 

Figure  3  illustrates  the  three-dimensional  model  of  the  first-order  solution  at 
point-moon  intercept  while  Figure  4  further  illustrates  aspects  of  the  three-dimensional 
model.  Figure  5  delineates  a  flowchart  of  tasks  to  be  accomplished  in  the  three- 
dimensional  model  of  the  first-order  solution  of  point-Moon  and  massless  spherical  Moon 
intercepts.  This  flowchart  summarizes  steps  used  to  construct  Appendix  B.  Symbols  and 
abbreviations  are  defined  in  the  Table  of  Symbols  and  Abbreviations.  The  hyperbolic 
trajectory  yielded  an  energy  requirement  of  63.2  kmVsec"  while  the  Apollo  mission 
required  a  minimum  energy  requirement  of  61.8  kmVsec".  This  results  in  a  2%  increase 
in  energy  by  using  the  hyperbolic  trajectory,  and  as  a  consequence,  will  increase  launch 
costs.  Considering  that  the  time  of  llight  is  reduced  by  more  than  50%,  the  2%  energy 
increase  penalty  seems  feasible  for  this  mission, 

2.  Three  Dimensional  Considerations 

Figure  6  illustrates  the  relationship  of  orbit  eccentricity  and  spacecraft  true 
anomaly  at  moon  intercept.  Figure  7  is  an  exploded  view  of  Figure  6  for  user 
interpretation  which  shows  eccentricities  between  1 .0  and  2.0.  These  diagrams  must  be 
referenced  by  the  user  to  insure  appropriate  inputs  for  solution  in  Appendbc  B.  Figure 
8  illustrates  a  cone  sweeping  out  a  circular  area  of  the  Earth.  The  Earth  surface  within 
the  cone  defines  a  zone  of  inaccessibility  for  Earth  launch  sites,  given  the  specified  true 
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Figure  3.  Three  dimensional  point  moon  model. 


Figure  4.  Three  dimensional  model  focusing  on  the  point 
Moon. 
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Figure  5.  Translunar  Trajectory  Program  Flowchart. 


anomaly  at  moon  intercept.  This  Tact  constrains  inputs  to  the  program  in  Appendix  B. 
If  the  launch  site  is  found  to  be  in  the  restricted  cone  region  a  different  trajectory  must 


be  used. 
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Figure  6.  Eccentricity  versus  maximum  and  minimum  true 
anomaly  at  Moon  intercept. 


Figure  7.  Exploded  view  of  Figure  6. 
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C.  SOLUTION  FOR  MASSLESS  SPHERICAL  MOON  INTERCEPT 

The  final  solution  in  Appendix  B,  is  the  massless  spherical  Moon  intercept.  An 
iterative  technique  is  employed  to  obtain  an  impact  on  the  lunar  surface.  This  technique 
is  illustrated  and  described  in  section  C.l.  Sclenographic  coordinates  are  utilized  to 
obtain  the  solution  for  the  massless  spherical  moon  intercept  and  are  described  in  section 
C.2. 

1.  Iterative  Technique  for  Lunar  Surface  Impact 

Figure  9  illustrates  the  method  to  ascertain  the  lunar  surface  impact  point 
iteratively.  The  iterative  process  is  defined  in  the  following  steps: 
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Determine  point-moon  intercept  coordinates  (R=pem>  fM>  and  Im  are  inputs). 


•  Select  trial  value  of  radial  distance  :  a  Pem'D. 

•  Calculate  the  corresponding  flight  time  from  perigee,  to,  and  separation  flight  time 
Do. 

•  Find  the  next  iteration  :  R,  =  Ro  -f-  (Do  -  D)cosAo. 

•  R„*D„-^R„D„  =  cosA.„  and  R„  =  R„.,-»-{D„.,-D)cosA,„.,  (bold  print  represents  vector 
quantities). 

•  Continue  until  n  =  10  or  |  D„.,-D  |  <  1  km. 


This  iterative  technique  is  employed  in  Appendix  B,  however,  solutions  are  calculated 


manually  by  the  program  user  without  a  terminal  loop. 


Figure  9.  Iterative  technique  model. 


1.3 


2.  Selenographic  Coordinates 

Impact  points  on  the  near-side  of  the  moon  are  expressed  in  selenographic 
coordinates,  for  which  the  zero  degree  longitude  and  latitude  points  at  the  Earth.  The 
North,  South,  East,  and  West  directions  indicated  on  Figure  13  for  latitude  and  longitude, 
are  for  an  observer  on  the  Earth.  The  lunar  librations  (Roy,  1988)  are  neglected  in  this 
definition. 

a.  Euler  angles 

Figures  10  and  11  express  Euler  angle  transformations  which  are 
incorporated  in  Appendix  B  to  relate  selenographic  coordinates  to  Earth-centered 
coordinate  systems.  The  first  Euler  angle,(d>),  is  zero  because  this  angle  is  not  required 
in  this  analysis.  The  second  Euler  angle  rotation,(©),  is  the  angle  i^  in  Appendix  B.  The 
third  Euler  angle  rotation,(’PM).  the  angle  <i>M  in  Appendix  B.  (Goldstein,  1951) 

b.  Lunar  Impact 

Figure  12  illustrates  a  selenographic  coordinate  solution  for  lunar  impact 
of  the  spacecraft.  Since  time  reversal  is  employed  in  Appendix  B,  the  Moon  moves 
backwards  in  its  orbit  for  the  translunar  trajectory.  The  sign  on  is  reversed  in 
Appendix  B  for  this  reason. 

The  selenographic  latitude  for  surface  impact  is  designated  6 
and  the  selenographic  longitude  for  surface  impact  is  designated  a.  Figure  13  illustrates 
these  latitudes,  longitudes  and  their  signs.  The  solution  for  a  has  a  factor  of  rr  added  to 
it  to  insure  that  the  positive  X-axis  points  to  the  Earth.  Figure  14  displays  a  locus  of 
impact  points  for  varying  time  offsets  (t„,f).  The  axes  units  are  in  radians.  Varying  t„f( 
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•  mCUNATION  OP  LUNAR 


Figure  11.  Third  Euler  angle  rotation. 
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Figure  12.  Lunar  centered  frame. 


SELENOGRAPHIC 

COORDINATES 


Figure  13.  Selenographic  coordinates. 
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is  equivalent  to  varying  the  input  value  of  the  lunar  sweep  angle  (<^m)-  By  providing  an 
offset  from  the  point-Moon  intercept  declination  (5^  off),  the  user  can  target  from  North 
to  South  (see  Figure  13).  This  factor  can  be  combined  with  (t^fj)  to  produce  an  arbitrary 
location  of  the  impact  point  on  the  lunar  near  side. 


Figure  14.  Locus  of  lunar  impact  points. 


D.  PATCHED  CONIC  APPROXIMATION 

The  method  of  calculation  for  the  massless  spherical  moon  intercept  model  can  also 
be  used  for  sphere-of-influence  calculations  if  the  radius  of  the  sphere-of-influence  is 
substituted  for  the  radius  of  the  Moon.  The  patched  conic  lunar  impact  point,  using  a 
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finite  lunar  mass  model,  can  then  be  approximated  by  an  iteration  process  identical  to  the 
process  used  in  Appendix  B  (Bate  ei  al.,  1971). 

From  the  final  patched  conic  iteration,  the  impact  point  for  the  Moon-centered 
conic  trajectory  can  be  calculated.  This  calculation  is  made  using  the  same  method  as 
Appendix  B. 
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III.  APPLICATIONS 


Fast  Earth-return  trajectories  have  numerous  applications,  such  as  for  requirements 
in  evacuating  a  lunar  base  in  an  emergency  or  the  establishment  of  a  lunar-based 
strategic  deterrent  force.  This  chapter  will  discuss  such  applications,  their  development, 
and  technical  feasibility. 

A.  STRATEGIC  PEACE  INI'I  IATIVE  CONCEPT 

The  ABM  treaty  prohibits  nuclear  weapons  in  space.  If  we  look  at  a  future  world, 
this  treaty  might  possibly  be  altered  or  not  exist.  The  Strategic  Peace  Initiative  is  a 
hypothetical  plan  to  place  nuclear  weapons  on  the  Moon  and  to  eliminate  them  from  the 
Earth.  The  reasoning  for  this  replacement  is  to  prolong  the  strike  capability  and  reaction 
time  that  now  exists,  thereby  increasing  stability. 

1.  Background  Information 

This  section  describes  a  hypothetical  Strategic  Peace  Initiative,  a  plan  to 
eliminate  the  risk  of  global  nuclear  war  and  simultaneously  promote  peaceful  uses  of 
space.  The  SPI  concept  could  be  a  step  toward  reducing  present  nuclear  confrontation 
as  well  as  "denuclearizing"  the  Earth.  (Wad.sworth,  1988) 

Rather  than  using  our  current  posture  of  SDI,  the  SPI  concept  can  encourage 
peaceful  exploitation  of  space  (Gray,  1985).  SPI  resources  could  provide  the  logistics 
support  for  a  new  era  of  nonmilitary  space  activities,  such  as  manned  space  exploration, 
servicing  geostationary  communication  satellites,  safe  disposal  of  radioactive  waste,  clean- 
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up  of  hazardous  space-debris,  and  eventually,  deflection  of  Earth-threatening  asteroids 
or  comets. 

Although  the  basing  of  deterrent  strike  forces  in  space,  particularly  on  the 
Moon,  has  already  been  proposed;  this  proposal  is  quite  unique.  The  novelty  of  an  SPI 
concept  is  that  by  a  new  treaty,  the  United  States  and  Soviet  Union  would  share  the 
"high  ground"  of  space.  As  with  both  the  INF  and  Start  treaties,  success  of  the  SPI  treaty 
would  depend  on  whether  it  serves  the  mutual  interests  of  parties  involved.  (Wadsworth, 
1988) 

The  free  world’s  offensive  deterrent  policies  have  prevented  a  global  nuclear 
confrontation  to  date.  Unfortunately,  this  policy  is  dangerously  unstable  because  of  the 
critically  short  decision  time  for  launching  a  counterstrike  to  a  perceived  nuclear  attack. 
SPI  would  eliminate  or  greatly  reduce  the  danger  of  accidental  nuclear  war  by  increasing 
decision  time  from  minutes,  in  the  case  of  an  SLBM  attack,  to  a  more  comfortable  two 
days. 

The  primary  justification  for  SPI  is  survival  in  a  nuclear  age,  consequently  the 
emphasis  in  this  proposal  is  on  deterrence.  A  secondary,  economically-compelling 
justification  is  that  much  of  the  SPI  investment  in  global  security  could  also  serve  the 
civilian  space  program.  They  both  require  similar  logistics  elements: 

•  A  heavy-lift  launch  vehicle 

•  A  space  tug  (space/lunar  transfer  vehicle) 

•  Earth/Moon  space  communications  and  tracking  networks 

•  A  manned  lunar  base 
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2.  SPI  Architecture 


The  minimum-energy  iransearth  strike  trajectory  has  a  4.7  day  flight  time. 
Doubling  the  injection  energy  yields  a  three  day  trajectory.  High-energy  strike  trajectories, 
taking  two  days  or  less,  still  permit  nuclear  strike  weapons  to  escape  the  Earth-Moon 
system  if  aborted.  Appendix  B  provides  such  a  time  of  flight  (approximately  1.9  days). 
Shortening  the  flight  time  from  4.7  days  to  2  days  results  in  only  a  two  percent  increase 
in  transearth  trajectory  energy,  which  is  equivalent  to  providing  the  additional  velocity 
increment  typical  of  an  ICBM.  Atmospheric  re-entry  velocity  does  not  increase 
significantly  between  4.7  day  and  2  day  trajectories,  which  simplifies  re-entry  vehicle 
design.  To  minimize  vulnerability  to  intercept,  the  trajectories  must  be  of  the  direct- 
ascent  type  without  lunar  parking  orbits,  in  contrast  to  the  more  economical  supply 
shuttle  trajectories  which  must  also  be  utilized.  (Wadsworth,  1988) 

The  SPI  communications  networks  and  surveillance  systems  would  preferably 
be  located  at  the  natural  synodic  satellite  points,  designated  LI  through  L5,  which  are 
stationary  with  respect  to  the  Earth-Moon  system  and  would  minimize  station-keeping 
fuel.  These  points  are  therefore  an  ideal  location  for  placing  these  communications 
satellite  systems.  (Farquar,  1970) 

The  two  opposing  lunar  ba.se  zones,  each  approximately  one  hundred 
kilometers  in  diameter  would  contain  the  strike  force  missile  silos  and  manned 
command/support  complexes.  These  zones  would  be  located  on  approximately  opposite 
sides  of  the  moon.  They  might  occupy  polar  sites  to  take  advantage  of  full-time  sunlight 
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for  solar  power.  Determining  optimum  zone  locations  requires  analysis  if  the  optimum 
strike  trajectory  missions  for  both  bases.  (Wadsworth,  1988) 

Similar  to  the  deployment  of  the  stages  of  SDI,  SPI  would  be  implemented 
in  stages  and  be  subject  to  treaty  compliance  verification  (possibly  a  United  Nations  role) 
at  each  stage.  Both  superpowers  would  deploy  the  following  five  stages: 


•  NON-NUCLEAR  DEPLOYMENT  :  Lunar  bases,  communications,  and 
surveillance  systems  deployed. 

•  FLIGHT  TESTS :  Test  strike  force  with  dummy  warheads  to  demonstrate  accuracy 
and  reliability. 

•  NUCLEAR  DEPLOYMENT  :  Equip  lunar  bases  with  nuclear  missiles  and  phase 
out  Earth-based  strategic  nuclear  forces. 

•  VERIFICATION  :  Surveillance  of  lunar  bases  and  Earth-Moon  space;  inspection 
and  monitoring  of  Earth-based  nuclear  device  manufacturing  and  test  facilities;  and 
inspection  of  translunar  cargos. 

.  STRATEGIC  REDUCTION  EVOLUTION  :  Bilateral  reduction  of  strategic 
nuclear  forces  to  minimum  required  for  asteroid  deflections;  increase  in  cooperative 
nonmilitary  space  ventures  by  utilization  of  lunar  bases. 


B.  SPI  TECHNICAL  FEASIBILHY 

Initial  deployment  of  SPI  relics  on  technology  that  has  or  will  soon  be  proven  by 
other  programs.  In  fact,  SPI  can  virtually  be  used  as  a  spin-off  of  developed  SDI 
technology.  Programs  such  as  TDRSS,  NASA  space  station,  OMV,  NASP,  HEDI,  and 
the  APOLLO  lunar  programs  are  examples  of  some  of  these  programs. 

Challenges  in  technology  do  exist.  These  challenges  include  a  precision  guided  re¬ 
entry  vehicle,  developing  "smart"  surveillance  techniques,  and  prevention  of  adverse 


physiological  effects  of  long-term  residency  in  low  gravity.  At  an  interval  of  deployment 
of  ten  to  twenty  years,  the  fully  deployed  SPI  deterrent  force  system  probably  would  cost 
comparably  to  the  proposed  SDI  system.  (Wadsworth,  1988) 

Under  the  SPI  treaty,  a  sovereign  lunar  base  zone  would  be  assigned  to  each 
superpower  for  basing  its  nuclear  strike  force.  These  zones  would  be  located  on 
approximately  opposite  sides  of  the  moon.  The  optimum  size  and  location  of  the  zones 
would  be  a  compromise  among  the  following  six  basic  requirements: 


•  Minimize  Surveillance  Cost:  Requires  minimizing  zone  size;  locate  bases  on  the 
Moon’s  visible  side  to  take  advantage  of  surveillance  redundancy  provided  by  Earth- 
based  sensors  to  back  up  the  spaced-based  sensors. 

•  Maximize  Nuclear  Survivability:  Requires  dispersal  of  hardened  missile  silos  within 
a  circular  area  of  at  least  one  hundred  kilometers  diameter. 

•  Minimize  Base  Vulnerability:  Requires  maximizing  separation  between  bases  to 
deter  pre-emptive  strikes  and  reduce  vulnerability  of  transearth  strike  forces  to 
interception. 

•  Minimize  Transportation  Costs:  Requires  selection  of  base  zone  locations  to 
minimize  fuel  cost  for  both  supply  shuttles  which  utilize  parking  orbits  and  strike 
forces  which  require  direct-ascent  trajectories.  (This  issue  will  be  difficult  to  resolve, 
because  if  both  strike  forces  were  to  lake  advantage  of  the  Moon’s  orbital  velocity, 
then  locations  near  the  lunar  poles  will  be  favored.  This  induces  a  strategic 
problem,  because  the  energ>'  requirements  for  striking  Northern-hemisphere  Earth 
targets  will  differ  between  the  polar  bases.). 

•  Avoid  Base  Overflight:  Requires  designing  supply  shuttle  parking  orbits  not  to 
overfly  opposing  base. 

•  Minimize  Communication  Costs  :  Requires  designing  base  communications  on  the 
near-side  of  the  Moon.  Far-side  Moon  bases  require  more  communication  system 
elements  because  its  orientation  with  the  Earth  does  not  meet  full-time  line-of-sight 
communications  link  requirements. 
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C.  SPACE  SUPPORT  AC  TIVITIES 


A  large  variety  of  space  support  programs  can  be  served  by  SPI.  For  example, 
lunar-based  astronomy  would  be  free  of  restriction  imposed  by  Earth’s  atmosphere.  A 
lunar  base  could  be  used  as  a  staging  platform  for  a  manned  expedition  to  Mars,  taking 
advantage  of  a  smaller  round  trip  propulsion  requirement  for  liftoff  from  the  Moon, 
rather  than  from  the  Earth.  The  servicing  of  geostationary  satellites  could  be 
accomplished  with  SPI  resources.  Transport  of  concentrated  nuclear  waste  to  a  safe 
dumping  site  in  the  Moon  would  he  unhampered  by  groundwater  problems  as  on  the 
Earth  (Rosen,  1981).  SPI  resources  could  be  used  for  debris  removal  from  low-Earth 
orbits.  This  is  an  increasingly  serious  hazard  to  space  missions.  Also,  safe  disposal  of 
space  junk  or  dead  satellites  in  higher  orbits  could  be  an  SPI  priority. 

Using  the  SPI  surveillance  network,  the  nuclear  arsenal  of  the  Moon  can  be  used 
to  deflect  Earth  threatening  asteroids  or  comets.  Collision  with  a  ten-kilometer  diameter 
asteroid  has  the  explosive  force  of  ten  million  one-megaton  nuclear  bombs.  With  an 
explosion  of  such  magnitude,  the  impact  dust  cloud  will  plunge  the  world  into  a  deep 
freeze  similar  to  "nuclear  winter"  but  much  more  devastating.  Fossil  records  indicate  that 
such  an  asteroid  could  account  for  extinction  of  the  dinosaurs  and  cause  for  the  Ice  Age 
(Malove,  1985). 

D.  SUMMARY 

The  Strategic  Peace  Initiative  Concept  is  a  revolutionary  plan  for  strategic  nuclear 
deterrence.  It  is  a  plan  which  provides  double  value;  a  solution  to  accidental  global 
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nuclear  war  by  increasing  reaction  lime  to  at  least  two  days,  and  promotion  of  peaceful 
uses  of  space  by  developing  permanently  manned  lunar  bases  and  their  infrastructures. 

The  uniqueness  of  SPI  is  that  it  requires  bilateral  treaty  agreement  between  the 
United  States  and  the  Soviet  Union.  Although  the  politics  of  such  a  treaty  requires 
extensive  thought,  verification  of  the  treaty  can  be  accomplished.  By  looking  to  the 
future  of  peace,  perhaps  we  can  be  at  peace  with  ourselves  and  with  future  generations. 
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rv.  CONCLUSIONS  AND  RECOMMENDATIONS 

This  final  chapter  summarizes  the  findings  if  this  thesis.  It  also  recommends 
research  for  further  work  concerning  this  topic. 

A.  FINDINGS 

This  thesis  utilizes  the  selection  of  an  efficient  and  operationally  meaningful  set  of 
inputs  and  terminal  parameters,  to  lead  to  a  transearth  or  translunar  trajectory  targeting 
solution,  which  rapidly  converges  in  an  iteration  process. 

B.  RECOMMENDATIONS 

This  thesis  lays  the  groundwork  for  an  accurate  patched  conic  method  of  designing 
fast  transearth  trajectories.  Recommendations  include  the  finalization  of  the  patched 
conic  approximation  and  the  iterated  impact  point  for  the  Moon-centered  conic 
trajectory.  A  development  of  the  SPI  concept  to  include  a  cost/mission  analysis  would 
also  be  very  useful. 
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APPENDIX  A.  PATCIIED-CONIC  BASED  ON  INERTIAL  FRAMES 


The  patched  conic  trajectory  can  be  formulated  in  terms  of  Earth-centered  and 
Moon-centered  coordinate  frames  which  arc  inertial,  being  attached  to  a  fictitious  Earth 
and  Moon  which  travel  at  constant  velocity  with  respect  to  the  barycentric  frame.  The 
Earth-centered  conic  trajectory  approximation  implies  the  lunar  gravitational  force  is 
ignored  so  that  the  Earth  is  not  accelerating  with  respect  to  the  barycentric  frame. 
Consequently,  after  the  instant  the  spacecraft  is  launched,  the  center  of  gravitational 
force,  the  Earth,  is  replaced  by  a  fictitious  Earth  which  is  imagined  as  moving  in  a  straight 
line  at  constant  velocity  in  the  barycentric  frame.  In  the  meantime,  the  massless  Moon 
continues  on  its  circular  orbit  in  the  barycentric  frame.  In  terms  of  the  translating.  Earth- 
centered  frame,  the  lunar  orbit  appears  as  a  compound  curve,  rather  than  a  circle.  This 
fact  must  be  taken  into  account  in  targeting  the  conic  trajectory  for  lunar  intercept  or 
crossing  of  the  lunar-sphere-of-inlluence.  The  same  considerations  apply  to  the  Moon- 
centered  conic  phase  where  the  Moon  is  approximated  by  a  fictitious  Moon  moving  in  a 
straight  line  at  constant  velocity  in  the  barycentric  frame.  For  greatest  accuracy,  the  line 
is  taken  tangent  to  the  Moon’s  circular  orbit  in  the  barycentric  frame  at  the  time  of 
intercept.  This  requires  an  iterative  .solution.  (Wadsworth,  1991) 
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AIM’ENDIX  n.  LUNAR  TRAJECIORY  PROGRAM 


INPUT  PAR/AIKIKIiS; 


S/L'  dc.cliriali:n  ;5l  cciiili  [.cinl 

5,,  -  '28.3 

180 

(Inptil  in  degrees} 

Oj..  =  0.493928178314395 

radisns 

Ime  cr/imely  cf  s/c  ot  jneer  inlcr.'0{'l 

Im  " 

180 

(inpul  in  degrees) 

f  =  2.8623399732707 

Indians 

inclinalion  of  Uic  moon  erf'it  plane 

(inpul  in  degrees : 

iw  -~*20 

^  180 

2G.5  >  i  >  10.!i) 

M 

=  0.349065850398866 

Indians 

ccoenlricily  of  Itc  Innsfucar  Injcclco 
(eccen  >  =  eceen  } 

ccccn  ■  1.01 

niiri 


gcoceilric  allilude  cf  lacnch  ^oinl 
(ollilDdc  in  kilomclcr.;) 


hy  -  1 


lonar  weep  angle  at  moon  inleroepl 

(input  in  degrees) 

^M  - 

♦m  = 

radians 

posigrade  or  retrograde  oitit  detemioaUoD 
posigrade  =  +1 
retrograde  =  -1 
“E_orrr  =  i 

3m_  •-  0  (  Dol  iterated  io  program  ) 


IN1CRHCDUTC  OUTPUl: 


6m  ■=  “in[sin[iM]-sin[^M]] 

6  m  =0.059426145347819 
(spacecraft  declioalioB  at  noon  ictercepl) 


“Ml 


cos[5mJ 


“Ml  ="0.986549223316463 


“M2 


[co»[iM]-8in[^Ml1 


“M2  =  0.163464460888851 


“M 


2*atan 


[^'“Ml  ] 
“M2 


UM  =  2.97739131857372 


2.96705972839036 

\ 
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INPUT  CONSMNTS: 


P£M  .=  384400 

km  :  carib-moofi  meat  disteice 

r£  :=  6378 

km  :  cerii)  meat  radios 

fM  =  1738 

km  :  mooo  meat  radios 

U£  :=  7.292 *10“* 

red /sec :  earth  angilar  rolalico  rele 

“M 


w 

Tio 


6)  jjj  =  2^1627109291352*10”*  rad/scc  : 

noon  n«en  oi1>ilal  nlalNtal  reU 


2  3 

k£  :=  398600  km  /sec  :  caitk  fravilalional  ooidacl 
2  3 

k]4  =  4903  km  /sec  :  mooi  iravitoUcnal  coi^sol 


OUTPUT  TlRUaeS  «d  fURNUUS; 

oibil  (or  Uc  Iraoslaiar  Trajectory) 

Rei  =  «■£♦•»£ 

R £1  =  6379*10^  km  ;  (s/c  redial  diiance  al  eartli  bioch) 
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RmI  =  PEM 

^M1  ~  3.844*10^  km  ;  (s/c  ndiai  didtoocc  at  moot  ifitereei^l) 

r  :=  (radial  distance  relic) 

^E1 

r  =  60^22887599937 

[r-[l>ccccn«cos[tM]]-l 

^  eccen 

fg  =  0.726288833059664 

^EM  =  *M"*E 

fEM  =  2.136051140211037 

«8"fEM  =  »«n[*EM] 
sigiifEM  =  0.84445214187218 

Sfem  •“  ^ 

“EMI  = 

“EMI  =-0.641452764362017 
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[c”[tEM]~»«"[^E]‘»'"[^M]] 


COs[  6|7  I'COSlSi 


NOTE:  for  Sfon 

if  sign  is  pcsilive  Uei  ose 
if  sign  is  negelive  ose  -1 


“EM2  =  0.767162532383019 


“EM  - 


2 


[»'«EM1  ] 
“EM2 


“EM  =  2J67186781944015 


NOTE: 

if  “EMI  =  i**®*  "EM 

is  »*«e_obt 

if  “EMI  "  fi'®*  “EM 

is  0 


“E  = “M’“EM 

a  £  =  0.710204536629705  ( s/c  litkl  aseeusioi  al  laonch ) 

«E1  =  ■^♦fM'^M 
«EI  =  -OJ38678825666912 

«E2  :=  «-fM~*M 
«E2  =  0.219826534971273 

fM_MAXl  =  «n[«E]*sin[aM]*«»[“EM]*c“[*E]‘“»[*M] 
fM_MAXl  =-0.53563101113312 

fM_MAX2  =  ij-fM_MAXl*] 
fM_>lAX2=0-8^5214187218 
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,  -  ,  ,  coa  fwl 

Ao  :=  asia  tanfftpl*u4 JmI"  f - r — j - rv 

1  I  ™  “J  [[co.[Je]-co.[8„])JJ 

^ot  =  1 J70796326794897  -  0.496321232060315i 

CM^nkl: 

oEmIN  = 

®^MIN  * -0«164201335016073  ♦  0.49632123206031S 


a£  nisi  be  eRal«r  Uen 


e«»nmi«  =0.974647175723671 


“«Kymin 


l  +  ecccn„i„]] 


( =  DiBiDtD  specific  CDeigy 
a(  knock ) 


“eigyinin  =  -1.023795999276086 


^E1  = 


cos[5e  *sin[lEMj 


=  0.421410460300705 


PE2  ■■ 


*K^em] 


/3£2  =  0.906870014913464 


Pe 


2*ataii 


1~<^E1 

^E2 


Iliehk  pakk  axinoUi  el  tanidi: 


^E  =  1.135796261808675 


i|  :=  co«[6e]**^^e] 

i|  =0.798478510707822 

«2  =  Jl-i? 

12  =  0.602023311789352 
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1  :=  2*atafl  — : — 
‘2 


0<  =  l<=ir) 


35 


Ce2  =  0.787491447196649 


Ce  :=  2* 


aUii| 


^~Cei 

Ce2 


Ce  =  0.906728172342672 


(srtep  A  koocb: 

ar>c  >  =  o) 

e 


-t-eccen'cos 

l-fcccen 


(Hiicee  4kUBce) 


Rp  =  5.57010652841956*10 


(l+eccen) 


(p«rigee  velocity) 


Vp  =  11.99319583652338 


Ve  = 


(l-eccen)*kE 

ll«ElJ 

Rp 

(s/cvekcily  el  leQidi) 


Ve  =  11.21107362995628 
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rifl  s  0.1040S2S84194718 

*M1  =«»[*m] 

fjU  :=«ii[rM] 

eccen*Ym  ‘IiMa 
"  l>[eecen*fnji] 

TM2  =  0.99457179716821 

y w  :=  2*aun  ■*  (KiU  paU  6a(le  a(  hood  iaiefcepi) 

“  L  ‘>'M2  . 

y,^  s  1.4665SS060009S38 


38 


^M1  =-0^145987662734 


^M2  ■= 


/JM2  =  0.799890488437212 


^M2 


(flitll  rsUi  azimoU  al  nooa  iDlarcepl) 


=  2Ji44799326S705S 


Xm  «  -3.785601002578928*10* 


(RccTiimilficoosnmTESor 

TBCKOONirPOWI-NOOtl 

mCSCEPT) 


Yji  - 


Ym  =  6J72482025241596*10^ 


Zm  •■=  Rmi  *«“[*m] 


Zm  =2^82996752157727  *10^ 


^  ^  (liM^nar  Ifajeclory  aifaineal- 

u  =  0.180439339283008  flf-^rieee) 
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Cm  •= 


Cm  =  3M2779312SS3709  (sneep  angle  at  oimd  kUrcepl) 


Xm'  = 


Vm 


cos 


•[[7M2~«“[®m]*‘2*=«»[Cm->m]]*«»[«M-»]] 


Pm], 


Xm*  =  -1^737958076395 


Ym* 


cos 


Vm 


*[[7M2"»“[*m]*‘2*«»{Cm"7m]]*»“[«M"Q]] 


[♦[‘1‘rMl  *«»[«M‘0]] 


Ym*  =  0J)48513058379955 

Zm’  ■=  V|4*i2*cos[Ciyi-yM] 

Zjf  =  -0X»S4S766973131 

„  .  [Xm*Xm’]*[Ym*Ym’]*[Zm*Zm’] 

Rm’  -  - - -  p - — - - 

Rmi 

(ndiei  aelocil;  e(  pciol-iDOOB  kleroepl) 

=  1.6<i|8023680688986 
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(ciml&r  oibU  period  (or  radios  R  ^ 

P 


Tp  =  4.137199S89247244*10^ 


*Ei  = 


l-coa[fE] 

»“[^e] 


f  El  =  0J7999703209726 

"  l-coaffi 

*Ml  -  -^ry- 

“"[*M 


(mi  »  7.115369722384208 


Ehl 


‘•J 

“«n-l 

I-*- eccen 

lUeccea 

Ejii  =  0.05361869041727 


Ei,2  •" 


1  + 


eccen  - 1 


1  ^^eccen 


1-  ^ 

'Vl^eccen 


^M1 


Imi 


E||2  -  1.103630544045282 
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r  -  -  r‘ 

ecccn ~  1  •  eccen~l  , 

E*hl  ■”  - - — *^E1  ~  • -  ^E1 

***^  IL  yl*eccen  J  <\|l+ ecccn 


£,1,1  »  OX>S3644386077804 


EsIi2  “  134171391602585 

*£1  = 

tEi  =0,089467283524084 

•Ml  •  ^L5 

2*«*(eccen-l) 

*M1  “  40.02754953820076 

I  £  :=  Tp  •!  £1  (eaHk  iMBCl  Hiekl  lime  reblwe 

lo  9erite<  raeeate) 

t£  s  370.144008646907 


g««»*E»hl]~^hl 

2*»*(ecceii  - 1)^ 
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*M  =  Tp*ti^, 


(^oial-nooD  inUrcepl  figtil  tint 
lebUve  (o  peiigee  paaaaie) 


tjjj  =  1.656019615080179*10* 

*EM  ■  ‘ E  (pfliil-BOflii  (reosiiBer  Irejedory 

lii(M  tine) 

tEM  »  1^231817499371  *10* 

Itns: 


Ro  -  382986J97662S  (uy  nditl  didfloce) 


OOlPVrS; 


Vo 


Rp 


^e] 


Vq  =  1^724659380983 


[Rp'tl+eccen)] 

“x  = - S- - 


Rx  -  0^29233197611653 


Iq  :=  2*aUn 


eccen  + 1  -Rjj  j 

|3' 

eccen-  1  +  Rj  j 

fo  =  2^1954051456415 
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Cq  =«+*0 

Co  =  3.042393390739423 
3o  :=  JBin[sm(i)*5iii[Co]] 

So  =  0.059657751937721 

AffMl  = 

At*  Ml  =-0.996881642139325 

tanfSM] 

A«M2  =lS(ir 

Aom2  *0078911289214964 

ff  »"Aomi  ] 

Abm  “  3062599237378809 

tin[8o] 

^•in(i)*Ua^Co]| 


Aa  1  -  -0.996857191524672 
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Aa2 


t«i[«o] 

tin(i) 


6a2  =  0^9219566431173 


^0  :=  2*stan 


Aoto  =  3.062289992046826 


“i  =  “M~^M 

a^  =  -0.08S20791880S089 

ao  :=  Aao^tti 


oo  =  2.977082073241737 


Yl  - 


>1  -  0.104293660037493 


eccen*y|  ‘sin 

*ol 

I] 

1  '•'eccen'cos 

t] 

>2  =  0.994S46S461S8S93 


y0  ■=  2*atan 


[i-n] 


Y2 


YQ  «  1.46631266S3401S7 
fil  :=  i2*cos[ao-»] 
fi  I  =  -0A00131267822715 


P2  ■= 


cos  6q 


P2-  0.799901S32303S08 


fiQ  :=  2*sUn 


^2 


Pq  =  2J14461S3046498S 


Xq  :=  Ro‘«»[*o]*«»[®o] 

(RCCTiimjlKCOOBI)lfUTfSOr 

Xo  =  -3.771433240671579-10*  yg.  OH  ITS  TUUKTORY) 

Yq  :=  Ro*a»[«ol*““[«o] 

Yo  =  6J60990806241751  *10* 

Zo  =  Ro*«“[*o] 

Zq  =  2^3455104046146*10* 
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‘MO  -  ’ 


siji[fo] 


f]^0  =  7.105421 1020S0813 


u 

(eccen-  1) 

•^MO 

(U  eccen) 

1- 

(eccen- 1) 

*fMO 

(1  >eccen) 

E|^  =  1.1017SS4S870149S 


®^8h3 


j_  (cccen-1) 
['vC  i  +eccen) 


!♦ 


(eccen- 1) 
■^(1  ♦ecccn) 


Ejjj  =  U38578545032142 


<M0 


,[ 


«=c”‘E5h3l-Eh3l 


[2’«*( eccen- 1)*^] 


tMO  =  39.8219787493874 

*0  =  '*’p*<MO 

to  =  1M751474124978*10’ 


*M0 
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ton  ;=  0  (lime  dfel  ■■>el) 


atq  ■=  tM”*o^*on 


Atq  =  850.4873830398719 


♦mi  ■=«»s[om  ]*««[«  m] 


^|yll  =-0.984807753012208 


♦M2  - 


»H‘m1 


♦M2  =  0.17364817766693 


:=  2**Un 


[■-♦mi] 

♦M2 


=  2.96705972839036 


♦mo  =  ♦M'WM'A^^O 


♦mo  ^  2.969323408665169 
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XmO  =  PEM'o»[^Mo] 


(RCCTANGUUfi  COORDINATES  or  THE 


XmO  =-3.787102316833942*10* 
Ymo  =  PEM'«»[iM]*»u»[^M0] 
Ymo  =  6.191940094679928*10* 

Zmo  •=pem**“«[«m]*“»[*mo] 

Zmo  =2JS3681886822944‘10* 


NOON  AT  TIKE  t  ) 

0 
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DMO  ■•=  ♦[Yo-YMof  ♦[Zo'ZMof] 

MSTiNCE  eerREN  m  spicccKin  and  the  lonar  center  at  twe  i 

0 

Dmo  =  1-738000000264939‘10^  km 

*  the  COAL  IS  TO  DRHE  THIS  OISTANCC  TO  EQUAL  THE  RADIUS 
or  THE  NOON  FOR  AN  WPACT  ON  TIE  SURFACE  OF  THE  HOON 

DMO'fM 

enw  s  2M9394446052611  ‘lO'^ 


Apo 


[[^MO'Xo]*Xo»[Ymo-Yo]«Yo»[Zmo-Zo]‘Zo] 

[Rq'^mo] 


Apo  =  0^1263979792074  J 

Kmart  ■=  Ko  +  DMO’ApO 


X  3J43986656315015M0’  ( Meter lo  .«  fcr  Mil  ol 

Ro> 
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X  =Xm-Xmo 


X  -  1S0.13142SS01429S 

Y  =Ym-Ymo 

Y  =  80541930S6166769 

Z  -Zm-ZmO 

Z  =  293.1486S334783i 

MPOTS: 

(DJURinaEliPWS) 

$  = 

♦m  =  ♦mo 

♦  M  »  2.96932340866S169 
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OVimS: 

(TflE  Mirmi  sisucms  is  IWNSeBEP  iis  lOlOHS:) 


Aj  Bj  Cl 
A2  D2  C2 
A3  B3  C3 


(iKni.ijWAmfiUNSun()N) 


Ai  :=  C0B[ifn4]»cos(4)-co5(^)*5in(*)*8in[»i^] 

A 1  =  -0.985198313432347 

Bi  :=  -•iii[*M]*coa(4)-cos(8)*siii(^)*co8[»M] 

B 1  =  -0.171418444777856 

C  j  ;=  »iB(8)*8iB(#) 

C,  »  0 

A2  :=  coi[+M]*«in(4)*a»(8)*coB(4)‘»ui[^M] 
A2  «  0.161080647624348 

82:=  -iinJi'j4]*8in(4)-*-cos{8)*'co5(4)*co8[^H4] 

82  =  -0.925783585143099 
C2  :=  '8iii(8)*co8(4) 

C  2  =  -0342020143325669 
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A3  :=  sin(tf)*»in[i'M] 


A3  =  0.058628S610S1S85 


B3  :=  8in(«)*cc»[»M] 


63=  H)336957668364338 


C3  :=  cm(0) 


C3  =  0.939692620785908 


A  := 


^2 

A3 


Bl  Cl 
B2  Cj 
B3  C3 


B  := 


X 

Y 

Z 


C  :=  A*B 


C  = 


-285.9729519600736 

-€21.7234494824876 

12€79404S1383525 


5.^ 


Co  =  -285.9729519600736 


Cl  =  -821.7234494824876 


Ci  =  12.87940451383525 


X||]i|  •=  Co 

Yhin  =  Cl 

Ziun  =  C2 


6  :=  aUui 


8  «  0i)148017S6974441 


SEiewGfUPflic  utmi  rxt  wace 

HPiCTWTHGIOON 
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